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Abstract 



A theoretical investigation is undertaken into the creep 
instability of two common shells of revolution, the infinitely 
long cylinder and the complete sphere. The linear and non- 
linear visco-elastic material idealizations are used in the 
various theoretical procedures. The critical buckling time 
is greatly influenced by the parameter a/h ratio, magnitude 
of initial imperfections and material selection- The spherical 
shell has a lifetime longer than the cylindrical shell having 
the same a/h ratio, initial imperfections and material con- 
stants, while the critical mode numbers are the same. 
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NOTATION 



w 


transverse displacement, positive outwards. 


u 


displacement in n direction. 

1 


V 


displacement in n direction. 
2 


M ,M ,M ,M 
1 2 21 12 


moments per unit length as defined in Figure 2. 


N ,N ,N ,N 
1 2 21 12 


forces per unit length as defined in Figure 2. 


e ,e 

11 22 21 12 


total strains. 


P ,P 

1 2 3 


exteimal loads per unit area as defined in Figure 2 


X ,x ,x 

1 2 3 


orthogonal coordinates as defined in Figure 1. 


2,. 8, 


transverse shear forces per unit length as de- 
fined in Figure 2. 


R ,R 
1 2 


Principle radii of curvature of a shell and 
radii of curvature of the coordinate curves 

n . 

1 -2 


a 


cylindrical shell mean radius. 


a* 


out-of~ roundness parameter defined by equa- 
tion (4.19) . 


h 


cylindrical shell mean thickness. 


S 


arc length. 


A,B 


defined in equations (2.5 b-c) . 


2 


distance measured along a normal from the 
middle surface of a shell. 


E 


Youngs Modulus of Elasticity. 


C,K 


material viscous coefficients. 


H 


creep exponent defined by equation (3.4). 


R 


spherical shell mean radius 


R* 


out-of- roundness parameter of spherical shell. 
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NOTATION (cont.) 



H 


spherical shell mean thickness. 


A 

n 


eimplification function (cylinder) . 


t 


time. 


t 

cr 


critical time. 


B 

mn 


amplification function (sphere) . 


C /D 
n ran 


maximxom initial imperfection; cylinder/sphere. 


n 


harmonic in circumferential (0) direction. 


m 


harmonic in circumpolar (([)) direction. 


1 2 


lines of principle curvature. 


K 

1 2 


change of curvature. 


K 

12 

* * 

K K 

1 2 


half total change of curvature in tvzist. 
principle curvatures = ^ 

1 2 

r 


e ,e 
1 2 


membrane s trains . 




shear strain. 


0 


poisson's ratio. 


a* 


effective stress defined in equation (3.5b) . 


a 


4 4 

— Ch or — CH. 


0,x 


circumferential and axial coordinates in 
cylindrical shell. 


4>,6 


circumpolar and circumferential coordinates in 
spherical shell. 


n 


4 4 

•j kh or ” kH. 




Ch^ ai^ 

9 9 • 


c 


Kh3 km 3 

9 or 9 . 
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NOTATION (contj 



X 

X 



A 



/s 

T 

C 



T 

C 



s 



T 

S 



s. . 

13 

( ) 

0 

( ),i 

( ‘ ) 

(“) 

( ) ’ 



-aaw/^ (time constant linear viscous cylinder) . 

-af^VC (time constant non-linear viscous 
cylinder) . 

-3Raw/25 (time constant linear viscous sphere) . 

-3R^>v/2^ (time constant non-linear viscous 
sphere) . 

Xt (non-dimensional time linear viscous cylinder) . 

Xt (non-dimensional time non-linear viscous 
cylinder) . 

^t (non-dimensional time linear viscous sphere) . 

At (non-dimensional time non-linear viscous 
sphere) . 

n(n + 1) 

defined by equation (5.12c). 

a * /a * 

0 

stress deviator, 
initial value of ( ) . 

li L 

9i 

1 ( )_ 

at 

dominant part of ( ) . 

perturbation part of ( ) . 
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1 . INTRODUCTION 



1. I Preface 

The problems of creep deformation and buckling of structural 
members while often addressed in aerospace and chemical engineering 
have, until recently, been aphonic design considerations in ocean 
engineering. The rationale being that the ocean, itself, offers an 
environment very unaccommodating to creep with its relatively low tem- 
peratures. Further, the general use of creep resistive materials for 
relatively short periods of loading made creep insignificant in the de- 
sign of ocean shell structures. Recently, the escalating ddmand man 
and his society have placed on the ocean in pursuit of its understand- 
ing and resources has resulted in a rapid advancement in ocean technology 
This new technology has generated an increased interest in the creep of 
shell structures. This manifestation of interest has been stimulated 
by many events, the greater use of creep prone materials, alxominum and 
titanium, the longer duration structures are required to withstand com- 
pelling hydrostatic pressures, and the need to predict volumetric changes 
due to creep of shells which have stringent volumetric design constraints 
such as the neutrally buoyant and compressible, free floating, current 
measuring SOFAR floats of the Woods Hole Oceanographic Institution. 

In the past, the ocean engineer was perplexed with a two part 
dilemma when asked to predict the creep performance of his design, first, 
find -an acceptable analytical model for creep deformation of shells and 
second, locate accurate material constants for low temperature creep. 
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The purpose of this paper is to solve the first dilemma of 



developing a simple analytic technique to predict creep deformation of 
two common shells of revolution, the cylinder and complete sphere, suit- 
able to oceanic applications and which will be useful in preliminary 
stages of design and may be adequate for final design. Finally, to dis- 
cover if creep buckling is a phenomenon that needs to be actively con- 
sidered in the design of long term ocean shell structures. 

1.2 Historical Review 

A literary search of the existing literature reveals that few 
published works have been transcribed which are directly applicable to 
conditions experienced in the ocean. Most papers deal with columns, 
plates or cylindrical shells loaded in axial compression. A number of 
investigators have attempted solutions for the creep behavior of thin 
cylindrical shells loaded \ander uniform external pressure. To the best 
knowledge of this author, no paper has been published dealing with the 
creep performance of a complete spherical shell loaded under external 
pressure. 

The author’s dealing with cylindrical shells use a variety 
of simplifications^ but do have an assemblage of assumptions common to 
all their analyses. These similar assumptions include: 

(a) thin shell assurt^tion, so that the radial stress compo- 
nent is negligible in comparison to the in-plane stress components 

(b) only the effect of secondary creep, according to Odquist's 
law (1937) are considered 

(c) buckling is not defined in the bifurcation sense, but 
rather as a state of progressive accelerated growth rate of the displacement 
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field, creating vinfavorable strain distribution and leading to full 
structural failure . 

(d) initial imperfections of a shell cause eccentricity of 
loading, thus serving as a lever arm for compressive in-plane forces, 
which in turn increases the bending moment, which in turn, increases the 
rate of curvature change 

(e) an assumed solution of ^ infinite trigonometric series, 
either a complete or partial Fourier series. 

Hoff, et al (1959) used a sandwich model, in which the outer 
skins support membrane and bending stresses and the core supports shear 
without deformation. The sandwich model avoids complications in the 
superposition of corrpressive and bending stresses due to the non- 
linearity of the creep law. Further, a quasi-elliptical initial irrper- 
fection was assumed. Samuelson (1970) developed a numerical solution 
for non-uniform external loads. Sankaranarayan (1969) generated a solu- 
tion for creep deformation of cylindrical shells under combined lateral 
and axial pressures for a material that adheres to the Tresca criterion. 
Berman, et al (1974) presented a numerical analysis for a cylindrical 
shell infinitely long, so that, a state of generalized plane strain exists. 

As did Hoff, Berman assumed that the initial in^erfection is quasi-elliptical. 

Penny and Marriott (1969) have carried out the only known, but 
rather special analysis of creep deformation of a spherical shell. They 
explored snap through of a spherical boss point loaded. 

All analyses errphasize the fact that initial imperfections 
play a most important role in creep buckling. That, buckling stems 
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from the accelerated growth of these initial inperfections in an otherwise 
iiniform initial displacement field, 

1. 3 Method of Analysis 

The conviction that creep buckling stems from the growth of 
small imperfections in the othervise uniform initial displacement field 
suggests the use of a perturbation method of solution as first developed 
by Grahams on and Goodier (1962) for dynamic plastic problems. The 
method assumes an axisymmetric dominant solution in which buckling does 
not occur and a perturbation solution in which the influence of bending 
is significantly greater than the in-plane forces thereby giving rise 
to inextensible buckling. 

Due to the axisymmetry of the dominant displacement field of 
the shell of revolution, only the transverse displacement field consists 
of a dominant and perturbation part, while the in-plane displacements 
have only a perturbation part so that 



w(n , 


n , 


t) = w(t) + 


w (n , n , t) 


(l.la) 


1 


2 




1 2 




u(n , 


n f 


t) = u* (n , 


n , t) 


(1.1b) 


1 


2 


1 


2 




v(n , 


n , 


t) = V* (n / 


n , t) 


(1.1c) 


1 


2 


1 


2 





respectively, with perturbation parts delineated by primes and dominant 
parts by bars. 

In turn, to remain consistent, the shear strain and twisting 

\ 

curvature change have only perturbation parts while the membrane strains 
and bending curvature changes have both dominant and perturbation parts. 
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similarly, the shear forces and twisting moments have just perturba- 
tion parts while normal forces and bending moments have both domincUit 
and perturbation parts 



N (n , 


n / 


t) = N (t) + 


N' (n , 


n , t) 


(1.2a) 


1 


1 


2 


1 


1 


1 


2 




N (n , 


n , 


t) = N (t) + 


N* 


(n , 


n , t) 


(1.2b) 


2 


1 


2 


2 


2 


1 


2 




N 


(n / 


n , 


t) = N* (n 


, n f t) 


= N (n , n , t) 


(1.2c) 


12 


1 


2 


12 1 


2 


21 1 2 




M (n , 


n , 


t) = M (t) + 


M' 


(n , 


n , t) 


(1.2d) 




1 


2 


1 


1 


1 


2 




M (n , 


n , 


t) = M (t) + 


M' 


(n , 


n , t) 


(1.2e) 


2 


1 


2 


2 


2 


1 


2 




M 


(n , 


n , 


rt 

11 


, n , t) 


11 

1 


(1.2f) 


12 


1 


2 


12 1 


2 


21 1 2 





The material comprising the shell is isotropic, incompressible 
and homogeneous. The initial elastic response is considered to be infini- 
tesimal. No initial shape of irrperfections is inferred. Further, the 
assemblage of assuitp>tions common to past analysis, (a) - (e) on page 11 
is assumed. 

The cinalysis will present non-linear equilibrium, strain- 
curvature relations and constitutive equations for linear and non- 
linear visco-elastic materials for an arbitrary shaped shell. The 
reduction of these relations will be carried out for the infinite cylinder 
and the complete sphere. 
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2. BASIC RELATIONS FOR THIN SHELLS 



2 . 1 Introducti on 

The succeeding relations were developed by VJalters and Jones 
(1972) . They were derived with the aid of differential geometry as 
discussed in Kraus (1967) and are valid for any thin arbitrary shell 
which has lines of principle curvature aligned with the parametric 
lines . 

2 . 2 Shell Coordinates - Differential Geometry 

The coordinates of point P on the middle surface of a shell, 

the reference surface, are defined in Figure 1. r is the position 

vector which is a function of parametric lines r , r . ^ is the unit 

1 2 

normal vector. The position vector can be expressed as 

r=x"x^+x^+xi^ (2.1) 

1 2 3 



where i, j, and k are the unit vectors in the x , x , x directions, 

1 2 3 

respe cti ve ly . 

As demonstrated in Kraus (1967) , the magnitude of the prin- 
ciple curvatures are the roots of the following quadratic equation 



- (EN-2MF-GL)k + (LN-M^) = 0 (2.2) 



where 



E = r, *r, (2.3a) 

1 1 



F = r , r, (2 . 3b) 

1 2 



(equal to zero when r is orthogonal to r ) 

1 2 
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2 



(2.3c) 



G = 



r, 

2 



• r 



H 



(EG-F^) 






(2.3d) 



/X — 

L = n*r, 

11 



(2.3e) 



— 

M = riT, 

12 



(2.3f) 



(equal to zero v/hen n 



v/ith 



cind n coincide with the lines of curvature) 
2 



N = n*r, (2.3g) 

22 



r, Xr, 






1 12 
’ll 



(2,4) 



Equations 5(a-d) are referred to as the first fundamental magnitudes 

while equations 5 (e-g) are called the second fundamental magnitudes. 

If n and n are orthogonal, then F equals zero, the 
1 2 

elemental area can be expressed as 



ds^ = A^d^n + B^d^n 
1 1 



where 



A = 



(E) 






and 



B = 



(G) 



(2.5a) 



(2.5b) 



(2.5c) 



2 . 3 Strain-Curvature Relations 

N A tensor formulation and the usual assumptions of shell 
theory, (eg. Kraus (1967)) have been employed by Walters and Jones (1972) 
to show that the strains in the thin shell may be written 
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e 



(2.6a) 



~ e + Sk 
11 1 1 



e = G + 2 k 
22 2 2 






and 



e + e = Y + 22 k 
12 21 0 12 






where t±ie membrane strains are 



u. A, V 
1 2 



G = 



w 

+ — + 



(w, )2 
1 



1 A AB R 2a2 



V, B, u 
2 1 



G = 



W 

+ + 



(w, )2 
2 



and 



2 B AB 2b2 

w, w, 

A.U, B.v, 1 2 

^0 ~ B^A , A B , AB 

2 1 



and the curvature changes are 



w, 



1 



A, 



K . i(.^i 

1 A R A , AB 'r B 

1 1 2 



w, B, 



1 , V 

= tt( 



1 



W, 



1 



2 B 'R B ^ AB ^R A ^ 

2 1 



and 



(2 












C 



2k = |( ^ 






12 BAR^ ^2 2 ABR^ g2 . 



c 



Note that y is the total shear strain while k is one-half the 
0 12 

twist curvature change 

\ 

The Love-Kirchhof f assumption has been invoked in the 
derivation of these equations. 



• 6b) 

• 6c) 

.7a) 

.7b) 

.7c) 

.7d) 

.7e) 

• 7f) 
total 
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2 . 4 Equilibrium Equations 



The typical shell differential element illustrated in Figure 2 



will have a rate of internal energy dissipation. 



W 



int 




N e + 


N Y 


+ 


M K 


2 2 


12 0 




1 1 


+ M K 


+ 2M 


fC 


Ids 



22 12 12 



( 2 . 8 ) 



where s is the total undefonned area of the middle surface of the 
arbitrary thin shell. Now, equating the rate of internal energy dis- 
sipation to the rate of external work (i.e., using the principle of 
virtual velocities) , a consistent set of equations of equilibrium is 
obtained. 



ABp - N B, + 
1 2 1 


(N A) , + (N B) , 

12 2 11 


+ N A, 
12 2 


(M B) , 

1 1 


M B, 

2 1 


(M A) , 
12 2 


M A, 
12 2 


T" " ■ " 

R 

1 


— . . ^ 
R 

1 


R 

1 


T 

R 

1 



(2.9a) 



and 



ABp 



N A, + 
1 2 


(N D) , + (N A) , 

12 1 2 2 


+ N B, 
12 1 


(M A) , 

2 2 


M A, 

1 2 


(M B) , 
12 1 


M B, 
12 1 

• 


R 

2 


— + 

R 

2 


CM 


*r 

R 

2 



0 



(2.9b) 



N AB N AB /N Bw, 



1 



ABp 

3 R 



R 



+ 



1 1 



fU Aw, 

2 2 



+ (N w, ), + (N w, ), + , 

12 2 1 12 1 2 L ^ J 1 



(M B) , 

1 1 
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Inertia terms are disregarded since creep is assumed to be a very pro- 
tracted deformation process 



3. CONSTITUTIVE RELATIONS 



3.1 Linear Visco-elastic Material 

The appropriate m.odel to simulate secondary creep during the 
application of a step function load is the Maxwell model illustrated 
in Figure 3a. The strain at all times is the sum of the elastic strain 
component, and viscous strain component. The stress in the viscous 
and elastic elements are always equal. Thus, from the equation of 
equilibrium the following expression for strain rate is found, 

(eg. , Drucker 1967) 



e 




(3.1) 



During application of the step function load, the motion of 
the dashpot is negligible. Hence, the initial strain is the elastic 
stretch of the spring. As time precedes, the model will continue to 
creep steadily at the rate of a/C as shown in Figure 3b. 

The relation for the linear Maxwell model can be written 
in tensor notation as, (eg., Drucker 1967) 



G . . 

13 



l+VA . (l-2v) . 3 ®ij 

~^ij T^i j ^ 2 " 



(3.2a) 
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where 



a 

KK 



S. . 
ID 



a 

11 



a. . 
ID 



+ (J + o 



22 33 



6 . . = kronecker delta 
ID 

and E, C are material elastic and viscous constants 

respectively. Since superposition of hydrostatic stresses is assumed 
to cause no further straining, equation (3.2a) reduces to 



e. . 

ID 




3!il 

2 C 



(3.2b) 



Inverting equation (3.2b) and performing the proper intergration reveals 
the linear viscous constitutive equations, after the initial elastic 
response, that is, disregarding the elastic behavior, for a thin shell as 



and 



where 



and 



N 




a(e 


.is ) 


1 




1 


2 2 


N 




a (e 


+ i e ) 


2 




2 


2 1 


N 




N 


a • 

= — Y 


12 




21 


3 0 


M 





C(K 


+ i K ) 


1 




1 


2 2 


M 





5(k 


+ - K ) 


2 




2 


2 1 


M 




-M 


= 


12 




21 


12 


a 


= 


|ch 






= 


Ch3 

9 





(3.3a) 

(3.3b) 

(3.3c) 

(3.3d) 

(3.3e) 

(3.3f) 
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3 . 2 Non-Linear Visco-Klastic Material 



Secondary creep for the non-linear visco-elastic material 
can be^ as in the linear case, represented by a Maxwell model where 



E 



(3.4) 



This model does not simulate primary creep and has a poor 
response to variable stress conditions. 

Assuming that the creep strains are dependent on loading 
history, the Prandtl-Reuss Incremental relations can be used for the 
non-elastic part 



where 



de. . 

13 



3 de* 

2 a* ij 



(3.5a) 



a* and e* are effective stress and strain based on Von 



Mses criterion 



where 

o* = — \(o -a )^ + (a -a )^ + (a -a )^ 

/2 L 11 22 22 33 33 11 

+ 6(t ^+T ^) 1 ^ (3.5b) 

12 23 31 J 



e* = — r (e -e )^ + (e -e )^ + (e -e )^ 

/J L 11 22 22 33 33 11 

+ 6(y (3.5c) 

12 23 31 J 



Norton's Law for secondary creep under a uniform constant load 
is expressed as 



de* = dt. 



(3.6) 
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Hence, one can substitute (3.6) into equation (3.5) and find 



e . . 
ID 



_ ^*M-1 

- S. 

2 K ij 



(3.7) 



which is Odquist’s Law for secondary creep, (eg., Odquist (1956)). 

By assuming that the effective stress remains constcint across 
the shell thickness, equation (3.7) can be inverted and intergrated to 
yield the non-linear constitutive relations while neglecting the 
elastic effects (Penney and Mariott 1971) 



N 

1 


* ( 1-M) • 1 • 

= Qo ^ ^ e ) 

11 2 22 


(3.8a) 


N 

2 


^ *(1-M) • 1 • , 

= (e + ~ s ) 

22 2 1 1 


(3.8b) 


N 

12 


- N - Y 

21 3 0 


(3.8c) 


M 

1 


* ( 1-M) • 1 • 

= CO (K + — K ) 

11 2 22 


(3.8d) 


M 

2 


^ *(1-M) , • 1 • , 

= CO (K + — K ) 

22 2 11 


(3.8e) 


M 


.. ^ *(1-M) • 

= -M = CO K 


(3.8f) 



12 21 12 

where 

4 

Q = 



Note, that if M equals one, the linear case, these constitutive 
relations reduce to the linear relations, equation (3.3). 
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4. INFINITELY LONG CYLINDRICAL SHELL 



4. 1 Introduction 

The objective of the following section is to reduce the 
general equations for an arbitrary shaped shell to the particular case 
of an infinitely long, thin cylinder. The perturbation solution for 
transverse deflection is assumed to be a Fourier cosine series 



w* 



CO 

A (t) cos n0 
n 

n=0 




(4.1) 



4.2 Basic Cylindrical Shell Equations 
4.2.1 Differential Geometry 

Defining x as n and 0 as r) , as illustrated in Figure 4, the 
1 2 

position vector equals 



— /\ /\ /X 

r = xi + asin0j + acos0k 



The principle curvature magnitudes are 
★ 

K = 0 

X 



and 

Kg = 1/a 



and the differential area is 



ds2 = dx^ + a^de^ 



where , 



A = 1 



and 



B = a 

Strain Curvature Relations 

Equations (2.7 a-f) reduce to 



(4.2a) 



(4.2b) 



(4.2c) 



(4. 2d) 
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u 



(4. 3a) 



e - 



'x + 



1 , . 
'^e ■ ? <''■6 + “ + 



^'0 'x '0 

Y = + V, + 

n a X a 



K = ~W, 

X XX 



= T ^^'0 - '^' 00 ^ 

a 



(4. 

(4. 

(4. 

(4. 



and 



2^0 = I - 2'^'x0^ 



(4< 



The time derivatives of (4.3a~f) are 



e = 



u, + w, w, 

X XX 



(4. 



"0 = a ^^'0 ■*■ i"''0'"'0^ 



(4, 



• 1 • • 1 • 

Y = - u, + V, + — (w, w, +w, w, ) 
na0 xa x00x 



(4, 



and 



K = “W, 

X XX 



'" d '~ '''■e ■ ”'eo’ 

a 



2k = — (v, - 2w, ) 

x0 a X 0x 



4.2.3 Equilibrixiiti Equations 

Equations (2.9 a-c) yield in-plane equations 



(4, 

(4 

(4 



ap + N . . + aN =0 
X x0 , 0 x,x 



(4 



and 



3b) 

3c) 

3d) 

3e) 

3f) 

4a) 
4b) 
4c) 
4d) 
. 4e) 

.4f) 

.5a) 
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0 



(4.5b) 



ap^ + aN +N +---M +M 

0 x0,x 0,0 a 0,0 x0,x 



and transverse equation 



ap - N- + a(N w, + N w, ) 
S 0 XXX x,xx 



I®e"'ee ”6,e"'e’ * 2”xe”'x6 



+ N w, +N w, + aM 

x0,x 0 X0,0 X x,xx 

M 

+ + 2M . - = 0 (4.5c) 

a x0,x0 

4.3 Buckling Solution 

For the axi symmetric loading of an infinitely long cylinder, 
the following assumptions are stipulated: 

(1) loading is axi symmetric, thus the dominant solution has 



no 0 dependence and p^ = p^ = 0 

(2) plane strain conditions exist, N 



N 



0 



X 2 

(3) transverse displacement, forces and moments do not vary 
significantly with the longitudinal, x, direction, hence both the dominant 
emd perturbation solutions are independent of x. 

The equilibrium equations simplify to dominant equation of 

n7 



P --1=0 

a 



and perturbation equation of 



(4.6a) 



'06 ‘‘‘ '■’'o,ee ° 



(4.6b) 



The dominant rates of membrane strains are 



e = 0 

X 



(4.7a) 
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i 




'I 




(4.7b) 



~ w 

^0 " a 

and 

and the rate of curvature change become 


(4.7b) 

(4.7c) 


K=0 Ka = 0 = ^ 

X 0 x0 





Assuming inextensible buckling the perturbation membrane strains 
are zero cind the curvature changes are simplified to 



K * =0 

X 


(4.7d) 


^'e = - -7 - “’'ee’ 

a 


(4.7e) 


K* = 0 

x0 

4.3.1 Linear Visco-Elastic Case 


(4.7f) 



Substituting the constitutive relations (3.3 a-f) into 
(4.6a) reveals a dominant transverse displacement rate of 



? 

• a p 

— 2 
w = 

a 


(4.8) 



The perturbation solution (4.6b) , with the assumption of 

inextensible buckling (e* = e'=e:' ^ = 0), becomes 

X 0 X0 



• 

- C/a(w-,99 + = 0 


(4.9) 



Substituting the assumed solution for w', (4.1 ) , into equation 

(4.9) yields a first order, ordinary differential equation of constant 
coefficients . 



XA + (n^-l)A = 0 

n n 

where 


(4.10a) 



- 26 - 



= amplification function 



Thus : 



where 



X = 



-aaw 



= C exp/-^^ 

^ ^ V(n2-1) / 



(4.10b) 



(4.11) 



Defining 



= an arbitrary constant dependent on the 
initial imperfections. 



= Xt 



one can plot A /C versus x to discover the critical mode (n . . 

^ n n c crit 

Figure 5 reveals that the elliptical mode, n equal to 2, is the critical 
mode. The growth of the critical amplification function is displayed 
in Figure 6. 

4.3.2 Non-Linear Visco-Elastic Case 

Combining the equilibrium equation (4.6 a~b) and the non- 
linear constitutive equations (3.8 a-f) generates a dominant transverse 

displacement rate of , , x 

2 ^.(M-1) 

a"^p a* 

w = ^ (4.12) 



The perturbation equation is transformed to 



- C/a(;-,QQ + w',QQQg) = 0 (4.13) 
Note, if creep exponent M equals one, equation (4.13) reduces to the 
linear equation (4.8). 

Substituting the assumed solution for w* , equation (4.1), into 
equation (4.13) yields a first order ordinary differential equation of 
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constant coef f icients . 



where 



XA + (n^+l)A = 0 

n n 



(4.14a) 



and 



hence 



= amplification function 




(4.14b) 



(4.15) 



v;here is an arbitraiy constant dependent on initial imperfections 
and 



T = Xt . 
c 



As in the linear case, the elliptical mode is critical. 

In the formulation of the analysis, it was assumed that the 
effective stress remained constant, however, as the perturbation solu- 
tion of transverse displacement grows, the effective stress increases. 

To overcome this problem one can divide the deformation process into a 
finite number of time segments in v;hich the effective stress is considered 
constant. This operation generates a piecewise approximation of the 
transverse deflection. 

Plane strain conditions allows the effective stress to be ex- 
pressed as 



a* = 0.866 ( 0 ^ + 0 * ) 

V u 

where 



(4.16) 



- 28 - 



(4.17a) 



which reduces to 



- 

h 



4 (1-M) • 

= ^^e'e) 



a. = 2 . 



(4.17b) 



Defining 



}i) = 



o* 



the perturbation solution becomes 



and 



. ^ / c , (M-1) . ^ ^ 

w* = C exp (“— \b ) cos2 0 

n 3 



w = -3 — w- 



(4.18a) 



(4.18b) 



where 



The effective stress becomes 



.866 



(■ 



l+6ip (H a*expl-~^ 



/t 



c , (M-in 






a* = (initial out of roundness pa- 

rameter << 1) . 

The initial effective stress is 



T = X t 

c 0 



X = initial value of X 



(4.19) 



p a 

.866 ■— (1 + 6 (J) a*) 



(4.20) 
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The duration of the time segment is dictated by the growth of 
effective stress. V/hen the effective stress increases by a given per- 
cent, i.e. ten percent (ip = 1.1) , a new time segment using the new value 
of a* is initiated. 

As an example, taking (a/h) equal to 10 and a* equal to .005, 
the amplification function's accelerated growth causes buckling at 
equal to 2.59, as graphically portrayed in Figure 7. Inspection of 
Figure 7 shows that the amplification function growth is very slow. It 
is only over the final 10 percent of the total lifetime that the ampli- 
fication function becomes substantial, indicating catastrophic buckling. 
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n 



5. COMPLETE SPHERICAL SHELL 



5 . 1 Introduction 

The succeeding chapter presents the creep response for 
the complete spherical shell. The perturbation solution for trans- 
verse deflection is assumed to be 



degree n and order m. 

5.2 Basic Spherical Shell Equations 
5.2.1 Differential Geometry 

Stipulating (J) as n and 6 as n ^ as exhibited in Figure 8, 

1 2 

defines the position vector equal to 



oo 



n 



w* 




(5.1) 



where P^ and P^ are Legendre and associated Legendre polynomials of 



r = R(sin<j) cos6 ^ + coscj) sin6 - sinc))'1<) , (5.2a) 



the principle curvature magnatudes equal to 



K 



= 1 
(f) R 



(5.2b) 



and 



* 

K 



_ 1 
6 R 



(5.2c) 



and the differential area equal to 



ds^ = R^d(})^ + R^sin^4)d6^ 



(5. 2d) 



where , 



A = R 
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" 9 ^ 





and 



B = R sin(|) 

Strain Curvature Relations 

Equations (2.7 a-f) simplify to membrcine strains 



■ -5 ("'♦ * ” * 



(5.3a) 



Eg = — (v,gCGC<^ + U COt(j) + W 



+ -^^5 — CSC ♦ , 



(5.3b) 



and 



0 R \ 0 9 



^^>^'0 2x 

+ — csc‘^<b 

R } 



(5.3c) 



and curvature changes 



and 



<J) ^2 ^ (J)(|) 

cscfj) , cotcb , , 

■"e = ~ “T- ^'^■"''<1,^ 



^ csccf) , . 1 , 

2k = (u, - W, ) + — (v, -V cot4> 

90 J^2 0 <!>0 J^2 9 



(5.3d) 

(5.3e) 



- w,g^csc<}) + 2w, gCSC<})COt(})) . 



(5.3f) 



The time derivatives of (5.3 a-f) are 



. I . . 1 

E. = t:(u,. + w + — w,,w,.) 
<}) R <1> R '<{) '(j) 



(5.4a) 



n 

• 1 • • * CSC^d * 

Eg = — (v,gCsc<^ + u cot<}) + v; + — — w,gW,g) (5.4b) 



* 1 / • • • csc^d) * 

Yq = -(csc<). u,g - v,^ + V cot<|, + — ^(w,^w,g 



+ w,^w,g)](5.4c) 
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and 



$ j^2 9 99 

• csc(j) • • , cot(J) ,• • . 

Kg = — (V, - W CSC<|)) + — (U-W ) 



2k = csc^ __ ^ i ^ ^ cot({) 

(j)0 ^2 9 99 j^2 9 



- w, ,_csc9 + 2w, csc9cot9) 

(pO 0 

5.2.3 Equilibrium Equations 

Equations (2.9 a-c) simplify into in-plane equations 



Rp9 - N0cot9+ N. . + N,cot9 

90,0 9,9 9 



+ ^ (Wj. i + M.cot9 - M^cot9 t M.n ^csc9) “ 0 
R 9,9 9 0 ^ ^ 



^ 2 COt <{. N^g + CSC < l . Ng^g 

+ i (CSC<t>Mg^g + M^g^^ + 2 M ^ g C O t <}) ) ^ 0 



and transverse equation 



+ N^(W,^COt<D + W,^^) + NgW,ggCSC^<j, 



- R(N^ + Ng) + 2csc<).N^gW,g^ + % , ' cj) 



+ csc^(t>N w, + csc<p(N w, + N qW, ) 



0 , 0 ' e 



90/9 9®'® 



+ + 2cot9M, ^ + csc^9M 

9/99 9/9 9 



0,00 



2m, Q -Cot(t)csc<t) + 2 m,_ ,.csc<J) - M. ,cot(j) 

( pO^tpO v 3 ^(p 



4 - Mg = 0 



(5.4d) 

(5.4e) 

(5.4f) 

(5.5a) 

(5.5b) 



(5.5c) 
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5 . 3 Buckling Solution 



In the analysis of the complete sphere, it is assumed the 

dominant solution has point symmetry, loaded in a uniform pressure 

field. As a consequence, there is no dominant 0 or (f) dependence, 

P0 = P4) = 0, M, = 0, M. = M-, cind N. = N,. The dominant transverse 
(J)0 (j) 0 0 9 

equilibrium equation simplifies to 

Rp - 2N^ = 0 (5.6a) 

2 9 

The in-plane perturbation equilibrium equations reduce to 



and 



cot(j)M* - cot(j)M* + cscc|)M* ^ 0 (5.6b) 

9,9 9 0 90,0 



csc(J)iyi* - + M._ . + 2cot({)M* ^ = 0 (5.6c) 

0,0 90 f 9 90 



The transverse perturbation equilibrium equation becomes 



N (w* , + cot({)w*, + csc^(J)w*, ) + M 



4) 



<t> 



00 (J) , ( 



+ 2cot(f)M* . “ + csc^([)M* + 2cot({>csc(j)M* . . 

9,9 (p 0,00 90,0 

+ 2csc(j)M* . - cot(|)M* , + M* = 0 (5.6d) 

<p0 , 90 0,9 0 

The dominant rates of membrane strains and curvatures 
changes are 



— — W 

^({) “ ^0 "* R 



e = 0 (5.7a-c) 

<p0 



and 



'^ 4 , = '^0 = ° 



K,. = 0 (5.7d-f) 

90 



Assuming inextensible buckling the perturbation rates of membrane strains 
are zero and the rates of curvature changes can be expressed as 
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(5.8a) 



K 




(5.8b) 



and 



• 2 

“ — (csc(j)cot(j)w’ , . - csc(j)w',,_) (5.8< 

(pd ^2 o 90 

If one examines the perturbation expressions derived for cur- 



(5.8c) 



vature changes v/hen inextensible buckling is assumed (equations 5.8a - 
5.8c), one discovers that they are identical to the curvature change 
expressions derived by Jones and Ahn (1974) . By coirparing the formu- 
lation of the transverse equilibrium equations of Walters and Jones and 

W 

Ahn, it is found that Walters and Jones assumed the — term to be neg- 
ligible in the curvature change expressions, while Jones and Ahn include 
it. The curvature change expressions are assumed to be a function of 
u, V, and w by Walters and Jones, while Jones and Ahn assume curvature 
change expressions are just a function of w. If, however, inextensible 
buckling is assumed, one can express the in-plane displacements, u and 
V as functions of the transverse displacement w only, and arrive at 
the expressions derived by Jones and Ahn. Thus, it is reasonable to 
assume that the transverse equilibrium equation of Walters and Jones 
must be modified by inclusion of the term 




(5.9) 



w 

which is generated by the — term. 

Hence, the modified transverse perturbation equation becomes 
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+ 2cot(J)M' 





The in-plane perturbation equations, (5.6 b-c) are rewritten 
in terms of w as 



if equations (5.10 b-c) are either equal to zero or nearly zero. The 
equations can be nearly zero since the transverse equation (5.10a) is 
considered the more dominant equation. 

5.3.1 Linear Visco-elastic Case 

The dominant transverse displacement rate is achieved by sub- 
stituting the linear constitutive equation (3.3a) into the equilibrium 
equation (5.6a) 




and 





To remain consistent, the preceding alteration is only valid 




(5.11) 



w 



3a 



Utilizing the following Legendre identities 




+ cot(l)Pn,^ + (csc^cj) + An)Pn = 0 (5.12a) 



and 



P^n, + no-hfhP^n. -h (m^csc^cf) + An)P^n = 0 (5.12b) 




v/he re 



An = n(n + 1) 
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The perturbation equation (5.10a) can be recast to 



where 



and 



$V^w' - (V^ + -I) (V^ + 2)w* = 0 (5.13a) 



( ) = ( ) + cot<j)( ) '(j, + csc^ ( )/QQ (5.13b) 



$ = 



3aaw 

25 



Substituting the assumed solution for the perturbation trans- 
verse displacement (5.1) uncouples the amplification function in the 
form of ordinary differential equations of constant coefficients 



4-AnB. + (t - Xn) (2 - Xn)B. = 0 (5.14a) 

On 2 On 

<i>XnB + (|- - An) (2 - Xn)B = 0 (5.14b) 

mn 2 mn 

Equation (5.14b) is identical to equation (5.14a) if m is understood to 
include the case m equals zero. 

Thus, providing p^ is constant 



B 



mn 



D exp 
mn 



^Ant 



(~ - An) (2 - An) 



(5.15) 



where is an arbitrary constcint dependent on initial imperfections. 

VThen is defined as a non-dimensional plot. Figure 9, can 
be constructed to exhibit that the critical mode is n equal to two. 

If the function 



w * 

D 

2m 



2 




ia=0 




(cos(j)) cosm0 



(5.16) 



is calculated for a number of points of various <f) and G values using 
computer program A of Appendix (A) , the two dimensional contour surface 
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I 







of Figure 10 is generated. On the contour map, positive values in- 



dicate an outward elevation of the mid-surface from the mean radius, and 
a negative value indicates an inward depression of the mid-surface from 
the mean radius, due to the perturbation solution. Since the elliptical 
mode has been proven to be the most critical mode, it will be henceforth 
assumed that the initial imperfections are of the elliptical shape. 

5.3.2 Non-linear Visco-elastic Case 

Substituting the constitutive equation (3.8a) into (5.6a) 
produces a dominant transverse displacement rate of 



2 

a^p a* 



w = 



3 ^ 



(5.17) 



In a method similar to the linear case (eqn. 5.12a), the 
perturbation transverse equilibrium equation is transformed to 



AV^w' - (V^ + + 2)w- = 0 (5.18a) 



with 



. 3 afb/ 

A - - j 



(5.18b) 



When the curvature change equations (5.4 d-f) are combined with consti- 
tutive equations (3.8 d-f) and then substituted into equation (5.10). 

Using the proper Legendre identities (5.13 a-b) , and assumed 
solution for w' (5.1), the uncoupled expressions for the amplification 
function become 



and 



AXnB + (t - Xn) (2 - Xn)B = 0 (5.19a) 

on 2 on 

AXnB + (t - Xn) (2 - Xn)B = 0 (5.19b) 

mn 2 mn 
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Thus, providing is constant 



Anx 



B = D exp , 
mn mn I ,3 



(5.20) 



where 



(~ - An) (2 - An) 



is an arbitrary constant dependent on 
initial imperfections, 



and 



T = At 
s 

The critical mode, as in the linear example, is n equal to two. 

As with the infinitely long cylinder, the assunption that the 
effective stress remains constant must be modified. Effective stress is 
assumed to be piecewise constant. Time is divided into a finite number 
of segments within which the effective stress is assumed to remain con- 
stant. When the effective stress increases by a given percentage, the 
segment is terminated, and a new segment with an increased value of 
effective stress is initiated. 

From equation (3.5b) the effective stress equals 



a* 




Oq) 2 + a 




(5.21a) 



where, from equations (1.2 a-f and 3.8 a-f) , the stresses are expressed as 



a, = + o^’ (5.21h) 

9 <P <P 



% = 



"o' 



(5.21c) 







(5.21d) 
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From equations (5.7 a~f) and 5.8 a-c) expressions (5.21b - 5.21d) can 
be expressed as 



P R 

0 = 0 = 

(1) 0 2H 



(5 



where H is the shell thickness of the sphere, 



rewritten in terms of w 



(5 



, 2 KH ^(1-M),3 , 

° A = “ T — ^ (-r w' + w' , 

<!> 3 p2 2 



” '(f) 

+ — - — COt({) 



w- , 



00 



csc^(J)) (5 



which can be further simplified using Legendre identities (5.12 a-b) 



KK (1~M) • 

a* (3w* 

3r2 






(5 



and in a similar manner 



KH (1-M) 3 • * • o 

o\ -a* i— w’ + w,^cot(}) + w' , qqCSc^(J) 






which is simplified to 



o' = a* (9w' + w, ) 

0 ■ 3^2 «•» 



*2 ‘ 



(5 



Finally the shear stress is written as 



, K ^(1-M). 

T ' = — O* Y 

4>e 3 ° ^12 



(5 



which, by substituting equation (5.8c), is 

,{1-M) 



T ’ 



2<H0* 



())0 



3R^^ 



(w ' , , csc(|)coti() - w'j.QCSccf)) (5 
(p (pU 



.22a) 

.22b) 

.22c) 

to 

.2 2d) 

5.22e) 

.22f) 

.22g) 

.22h) 
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J 




m 



since the perturbated strains depend only on curvature change 
rates (5.8 a~e) , which are all expressed as functions of w* , it is 
postulated that the spatial point on the sphere which exhibits the maxi- 
mum perturbated tremsverse deflection rate will be the point of maximum 
strains, thus stresses. To locate the point of maximum curvature rate (w*) 
one can utilize a precedure to study maxima and minima of functions of 

several variables outlined by Hilderbrand (1962). 

• • 

(a) Evaluate w*,, and w*,^ and then set equal to zero. 

9 0 

(b) Evaluate w',^^ , v> , and w’ , 



minimum exists. 



, and w* > 0 , ^a relative 



(ii) if w*,,w*_ > w*^ , and w* . < 0 , a relative 

99 00 90 99 

maximum exists. 

From equation (5.1) the critical value of w* (n = 2) is 

w* = B (t) fp (cos^) + cos0P^ (coscj)) 

2 m ^ o o 



+ cos20P^ (cos(j))^ 



with 



P (cos(f>) = •— Ocos^cJ) + 1) 

2 2 

dP2 (cos4>) 

P^ (coscf)) = sind) — — = Ssinicosd) 

2 d(cos(|)) 



P^ (coscj)) - sin^cj) 
2 



Thus , 



d ^?2 (cos(j)) 
d (cos(|)) 



= Ssin'^cfi 



(5.23a) 

(5.23b) 

(5.23c) 

(5.23d) 



• O'} 

w*,(j) = -3cos(j)sin(j) + 3cos0 (cos"^(j) - sin"^(j)) 

+ 6cos20sin(j)cos({) 



(5.23e) 
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(5.23f) 



• o 

v;* . = -3sin9sin<j)cos(j) - 6sin20sin"^cj) 

0 

Both of these expressions (5.23 e-f) are zero when <}) and 0 equal 
7r/2- 

Evaluating the second derivatives 



• 9 9 

w*, = 3(sin (}) - cos^cj)) - 12cos0sin(J)cos<J) 



+ 6cos20 (cos^(|) " sin^cj)) 



(5.23g) 



= -3cos0sin({)cos<J) - 12cos20sin^(!) 



(5.23h) 



and 



w\ = +3sin0 (cos^<J) - sin^c^) - 12sin2 0coscj)sin(j) (5.23i) 
90 



at 9 “ an<^ 0 = 7 t/ 2 one finds 



"■'oe = 



(5.24a) 

(5.24b) 






(5.24c) 



Thus, and a relative minimum exists at 9 

99 9“ 

equal to tt/2 and 0 equal to tt/2. Noting, however, that the deflections 

rates herein are of a negative sense (B is always negative) , the point 

2m 

9 = 'JT/2, 0 = tt/2 is actually a relative maximum. Further, this conclusion 
is in agreement with the maximum of the deflection field illustrated in 
Figure 10, which was calculated numerically. 

Substituting values (5.24 a-c) into equations (5.22 d, f, h) 

yields 

V 2 ^2 2m 






(5.25a) 
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(5.25b) 



, 7 KH ^(1-M) • 

o' = - — 0 * ' 'b (t) 
0 2 



2m 



<p0 



= 2 



KH ^(1-M) 






B (t) 
2m 



(5.25c) 



From equation (5.15) 



2 T X 

B' (t) = — -p exp 
2m 3 t 3 






the perturbated stresses simplify to 



(5.25d) 



^ p T 

. 3 , R, o _ . 2 , S , 

“ ♦ ' - I (j)"R- — exp(— ) 



(5.25e) 



o'e = a 



"'♦e " 3 (£)r* 



where 

R* = 

Hence equation (5.21a) becomes 



(5.25f) 



(5.25g) 



maximum deviation from radius of imperfections 

radius 



<1* - — H80' 

L ' 



)2 + ( 0 ^ + 0 ' )2 + (a^ - 7a' 



+ 6(-40j 



4> 



2]’’ (5, 



26a) 



simplified to 



0* = ^ r 

JT L 



, 20^2 _ 120 ^ 0 ’ , + 2000 

^ L 4> <«> 4> 






(5.26b) 



Substituting expressions (5.22a and 5.25e) , the effective 
stress equals 



0 * 



1 PR 






H 



1 9 ,Rv^. , s, 

I - 7 W* 

^ -ih 



+ 112.5(^)2R*2exp2(-^) 

H 3 _ 

The initial effective stress (x^ = 0) equals 



(5.26c) 
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1 



a * 
0 



(5.26d) 



Defining 



and 



1 PR r 1 

/2 



II R* + 112.5 (|)2r*2 

Z ri fi 




0 



T = A t 
S 0 



A = initial value of A (eqn. 5.18b) 
0 

the effective stress can be recast as 



o* = 



H 



T Ip 

1 9 ,R. ^,M-1 

exp — 



M-1 



+ 112(|)2R*2,|;(''^"l)^exp2 



M-1 



(5.27a) 

(5.27b) 



(5.27c) 



As an example, with R/H, R*, M, ip equal to 10, .005, 3, 

0 

1.1 respectively, the growth of the amplification function is as shown 
in Figure 11. Note that the amplification function stays relatively small 
until the last ten percent of the lifetime of the sphere. As expected, 
the time segment duration decreases as time increases. 

It should be further noted that the in~plane equilibrium equations 
(5.10 b & c) are equal to zero at the point of maximum stress and trans- 
verse deflection (c|) = tt/2, 0 = tt/2) . 



- 44 - 



6. DISCUSSION 



6. 1 Infinitely Long Cylinder 

6.1.1 Comparison to Past Results 

The most convenient basis of comparison of the analysis herein 
with past investigations is the predicted critical time of creep buckling. 
The time behavior of the initial irrperfection growth is such that it will 
become large and grow rapidly only for time near to the critical time. As 
cin aid for the reader. Appendix B has been provided to transform the nota- 
tion of past analyses into the notation used herein. 

Hoff et al (1959) predict a critical time of 



It is at this critical time at which the eimplifi cation of the 
initial imperfections becomes indefinitely large. The herein results have 
been compared to Hoff in three ways by varying the values of a/h, a*, and 
M in Figures 12 , 13, 14 respectively. The comparison reveals favorable 
agreement for smaller a/h ratios (a/h < 25) with a* and M equal to .005 
and 3 respectively; and small imperfections (a* < .003) with a/h equal to 
10 and M equal to 3. For large a/h ratios, the Hoff estimate of buckling 
time is approximately cin order of magnitude greater th£in the presented 
cinalysis. The Hoff solution is less sensitive to larger, gross initial 
imperfections, predicting a collapse time of five times greater than the 

V 

herein analysis when gross imperfections are assumed. 

The plot of varying creep exponent, M, (Figure 14) shows excellent 
agreement for M equal to one and seven and an estimate five times less than 
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cr 




( 6 . 1 ) 
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Hoff for M equal to three and five. In Figure 14 , a* and a/h equal 
.005 and 10 respectively. 

The actual time history of the two amplification curves are 
nearly identical, as exhibited in Figure 15. The amplification of im- 
perfections remains very small until the last ten percent of the life- 
time of the cylinder. 

The discrepancy in predicted critical time of the two cinalyses 
do not seem too serious in the light of the fact that creep buckling 
times are of the order of 10^ hours, so that excellent agreement is ob- 
served for the first 10^ hours. 

Berman et al (1974) presented the particular case of stainless 
steel 304, v;hich has the material properties at 1050^F listed in Ap- 
pendix C. One finds a great overestimation of Berman when compared to 
both the analyses of Hoff and herein. Berman predicts a buckling time 
for cylinder of a/h equal to 10 and a* equal to .05 of 1.743 x 10^ ^p^ 
while Hoff predicts 3.6852 x lO^^p^ 5-6789 ^ herein investigation 

estimates 1.826 x lO^^p^ 5-6789 ^ Both the Hoff and herein analysis 
predict a buckling time of approximately four orders of magnitude sooner 
than Berman. The author feels that the calculated time of Berman is in 
possible error, since he uses an out-of-roundness ratio of .05, a very 
gross out-of-roundness for any realistic engineering design. Further 
Berman fails to update the effective stress with time. The herein 
analysis does predict a buckling time of 8.3 x lO^^p^ 5-6789 
.005, (an order of magnitude less) and the effective stress is not up- 
dated with increasing time. 



The numerical analysis of Samuelson (1970) deals with the creep 
response of finite length cylinders subject to both non-uniform and 
uniform loads. Thus, higher harmonics come into an influencial role in 
the creep process. Further the transverse deflection is dependent on 
the longitudinal, x, coordinate 



w’ 



oo 

■E 

n=0 



ttx 

A (t) cosnGsin 

n L 



( 6 . 2 ) 



Samuelson predicts a buckling time of approximately 100 

seconds for a cylinder with a/h, L/a equal to 27 and 2 respectively. 

A similar buckling time is found for the same cylinder by the herein 

K 

analysis if M equals 3, a nominal a* .002, is assumed and — is 
assumed to be of the order 10®. Thus, it is possible to state that the 
herein analysis predicts a buckling time of the same order of magnitude 
as Samuelson. 

6.1.2 Elastic Response 

The elastic response of the cylinder can be obtained by sub- 
stituting the elastic part of the constitutive relation (3.2a) into the 
transverse equilibrium equation (4.6a) to yield 



w = 
e 



a^p (1-0^) 
g 

Eh 



(6.3a) 



The elastic transverse deflection from Roark (1965) is 
a^p 

w = (1 - -) 

e Eh 2 

Roark 



(6.3b) 



If Poisson's ratio, u, equals .3, the ratio of the herein 
elastic displacement to that of Roark is 1.07. 
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6.1.3 Parametric Study 



Due to the dependence of the time history of the transverse 
deflection on a multitude of independent variables, radius to thickness 
ratio (a/h) , out- of- roundness ratio (a*) , material constants (M, k) and 
loading (p^) / ^ indefinitely large number of parametric curves are 
required to adequately explore the creep response of cylinders under 
external pressure. This task is greatly simplified by the computing 
of transverse displacement (4.18a) using computer program B of Ap- 
pendix A. 



Figures 16, 17, and 18 are plots of non-dimensional critical 

't jp rn 
CITB 



time, [ - t .- t - versus a/h for creep exponent M equal to 3, 5, 7 re- 
spectively. In each case as a/h ratio for a given a* decreases, the 
non-dimensional critical time drastically increases. Further, as a* 
increases for a given a/h ratio, the non-dimensional critical time de- 
creases. Of the two parameters, a/h and a*, a/h is the more dominant 
variable. That is, critical time is more sensitive to a/h than a*. The 
same conclusion v/as reached by Hoff (1959) . 

Figure 19 illustrates the effect of veirying the initiation 

point of a new time segment of constant effective stress. If for example, 

t P ^ 
cr^ 

M, a/h, a* equal 3, 10, .005 respectively, tlie value equals 

3.5 X 10 ^ for time segments which are initiated when the effective 
stress (4.19) increases by fifty percent (i|) = 1.5) compared to a value 
of 2.2 X 10 ^ v/hen the time segments are initiated when the effective 
stress increases by five percent (tj) = 1.05), a ratio of 1.6 difference. 

From a practical ocean engineering point of view, a comparison 
of the two most common creep prone materials, aluminum and titanium, used 
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in the ocean is in order. The alloys used in the comparison are Aluminum 
6061-T6 and Titanium 6Al-2cb~lTa- . 08Mo. The room temperature creep 
constants of these materials are listed in Appendix C. The titanixom con- 
stants where experimentally found by Chu (1972) , while the aluminum con- 
stants by the author. 

The first comparison is non-dimensional critical time verses 
the initial elastic hoop stress, expressed as a percentage of yield stress. 
The yield stresses of Titanium 621/. 08 and Aluminum 6061-T6 are 120 and 
40 ksi respectively. The two materials are compared with two a/h ratios, 

10 and 50; and two a* ratios, .01, and .001. A number of conclusions 
can be drawn. The lifetime increases for both materials, approximately 
five times for a/h equal to 10 and forty-five times for a/h equal to 50, 
with an order of magnitude decrease in initial out-of-roundness ratio 
(a*). Hence, it appears to be more important that cylinders of large 
a/h ratios to be constructed to closer tolerances than cylinders of a 
smaller a/h ratio. Titanium is the better performer at low stress levels, 
but is surpassed by aluminum at higher stress levels. Imperfections play 
a more important role at higher stress levels. Finally, for a stress 
level less than forty percent of yield, there is a dramatic increase in 
lifetime . 

The influence the varying a/h at a given stress level, sixty 
percent of yield, has on the non-dimensional critical time is portrayed 
in Figure 21. Again, a crossover occurs. For small a/h ratios (a/h < 20) 
titanium is the optimum material, while for larger a/h ratios aluminum 
exceeds titanium. 
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TaJcing another tack, if a* is varyed for a consteint a/h (10) 
and a constant stress level, fifty percent of yield, Figure 22 is gen- 
erated, a crossover of performance is observed with titanium being 
surpassed by aluminum for grosser tolerances (a* > .007) . 

The most pertinent collation is the comparison of the two 
materials at equal pressure levels. As an example, with p , a*, a/h 
equal to 2000 psi, .001 and 10 respectively titanium out performs 
aluminum by four orders of magnitude, (Figure 2 3) . Titanium, with its 
superior strength, is able to operate at a much lower percentage of 
yield stress. 

Complete Sphere 
6.2.1 Elastic Response 

By substituting the elastic part of the constitutive equation 
(3.2a) into the transverse equilibrium equation (5.6a), an expression 
for the elastic response is developed 



w = 
e 



a^p^d - u^) 
3Eh 



(6.4a) 



The elastic deflection given by Roark (1965) is 



w 



a^p (1 - 0) 
g 

2Eh 



(6.4b) 



Roark 

A ratio of the elastic deflection to that of Roark is .86 of u is 
assumed to be .3. 

6.2.2 Parametric Study 

As with the infinitely long cylinder, the dependence the 
time history of the transverse deflection on a large number of indepen- 
dent varicibles requires tlie generation of a large assemblage of para- 
metric curves to adequately investigate creep response of the sphere. 
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Again, the task is simplified by numerically computing equation (5.1) 
using computer program C of Appendix A. 

Figures 24, 25, and 26 are graphs of non-dimensional time 
versus R/H for various creep exponents, (M equal to 3, 5, 7). 
The non-dimensional critical time decreases rapidly with both increas- 
ing R/H and R*, being most sensitive to R/H. 

The non-dimensional critical time of the sphere is fairly 
non-sensitive to the duration of the time segments, as illustrated in 
Figure 27. If the updating of effective stress (5 . 26a) is carried out 

when effective stress is increased by fifty percent (ij; = 1.5) the 
t P ^ 

cr^ 2 — R 

value of — is 1.6 X 10 If, however, the effective is updated 

t p 2 

C27 g 

when it has increased by five percent, = 1.05) — equals 1.05 

K 

X 10 an increase of fifty percent. 

When the non-dimensional critical time is calculated for 
various initial elastic loadings, expressed as percentage of yield 
stress, for the two alloys. Aluminum 6061-T6 and Titanium 621/. 08, 
Figure 28 is manifested. Initial imperfections are quite influencial 
on the critical time, causing a two order of magnitude increase for 
a/h equal to 10 and six orders of magnitude increase for R/H equal 
to 50 with an order of magnitude decrease in initial imperfections. 
Thus, revealing that for spheres of large R/H ratios, close tolerances 
are a prerequisite. Except for the case of high R/H ratio and gross im- 
perfections, titanium is the superior material, producing a longer 
lifetime. Rapid increase in lifetime is achieved if the stress level 
of the sphere is less than .sixty percent of the yield stress. 
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When R/H is varied for a given stress level, sixty percent of 
yield, titanium, regardless of the R/H ratio, out-perforras aluminum. 

This superiority is even more obvious, as demonstrated in Figure 29, 
for substantial R/H ratios (R/H > 30) . 

If for a constant R/H ratio and a loading of fifty percent 
of yield, the magnitude of initial imperfections is allowed to veer, 

Figure 30 is developed. Again, titanium is always the superior per- 
former, especially with small initial imperfections. That is, one 
has much more to gain by machining a titanium sphere to exacting 
tolerances than an aluminum sphere. 

At the same pressure level, the titanium is vastly superior 
than the aluminum, having a lifetim.e of approximately 7 orders of 
magnitude longer with p^ equaling 2000 psi and R* equaling .001. The 
obvious reason being that a titanium sphere will operate at one-fourth 
the percentage of its yield stress at a given ambient pressure. 

6 . 3 Comparison of the Creep Performance of the Sphere and Cylinder 

If the critical time of a cylinder and a sphere of the same 
material at equal pressure levels are compared by plotting the ratio of 
critical time of the sphere to the critical time of the cylinder for 
various a/h ratios, initial imperfections, cind creep exponents, (Figures 31, 
32, and 33) it is discovered that regardless of the initial conditions, 
the cylinder has tlie inferior creep resistant properties. This point is 
especially manifested in the case of titanium, M equals 7. The oscil- 
lations observed in Figures 31, 32, and 33 result from the variation of 
the ratio of the effective stresses with various a/li and a"*" combinations. 
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The sphere is more sensitive to variations in initial im- 
perfections, which is to be expected since the cylinder is a develop- 
able surface and is less constrained than the sphere. 

The sphere has excellent creep characteristics for stress 
levels less than sixty percent of yield stress. The cylinder, on the 
other hand, the stress level must be less than forty percent of yield 
stress for favorable creep properties. 

Finally titanium, because of its high yield stress is in 
every instance the best material to be used in ocean engineering ap- 
plications of both the sphere and the cylinder. 

^ Design Criterion 

The extensive analysis present suggests an assemblage of 
simple design guidelines to assist the engineer in developing a cost- 
effective, creep resistant design. In most ocean engineering shell 
designs, the major constraint in design is considered to be the minimi- 
zation of the weight to displacement ratio of the structure, while 
maintaining a maximum stress level of approximately fifty to sixty per- 
cent of yield. It is this constraint which dictates the use of light- 
weight, but creep-prone materials in the first place. Fortunately, 
this criterion is satisfied if titanium is used rather than aluminum in 
pressure vessels, even though titcinium is 1.5 times as dense as aluminum. 
The reason being, the aluminum shell would have to be four times as 
thick to achieve a stress field of equal percentages of yield stress. 
Hence, titanium will always give a longer lifetime. 

Wliile the sphere furnishes a vastly superior creep response 
than the cylinder, it is more expensive to construct and machine. Thus, 
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the engineer must weigh the merits of the costs of materials, con- 
struction and machining of both an appropriate sphere or cylinder for 
his particular situation (pressure level and time constraints) before 
he develops a final design. 

Regardless, the shape of the amplification curves for both 
the shell shapes suggest a creep design limit. For approxim.ately the 
first seventy percent of the lifetime imperfections remain small and 
grow in a linear manner. Thus, an arbitrary design limit can be estab- 
lished which states that as long as the imperfection growth becomes 
linear, the design will not creep buckle. The in-service monitoring 
of the creep deformation process is easily accomplished v/ith a three 
element strain gage rossette placed at the point of greatest imperfec- 
tion. When the effective stress approaches the magnitude of the yield 
stress, the active use of the structure should be terminated, 

6.5 Example Calculations 
6.5.1 DSRV ALVIN 

The research submersible, ALVIN, owned by the Office of 
Naval Research and operated by Woods Hole Ocecinographic Institution, 
has a spherical titanium 621/. 08 alloy pressure hull. The 40.005 inch 
radius hull has a a/h ratio of 20 and a maximum out-of-roundness ratio 
(a*) of 0.002. The maximum operating depth is 12,000 feet (5280 psi) 
and a collapse depth of 18,000 feet (7920 psi) . 

The secondary creep law for Ti 6Al-2cb-lTa-0. 8Mo alloy at 
room temperature has been empirically measured by Chu (1972) (Appendix C) . 
While the measurements were done in tension, it is believed that the 
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compressive creep behavior is siitu.lar to that of tension. The law is 
expressed as 

= Ke (6.2a) 



where 

y 

K = 8.3941934 x 10^° hr. (6.2b) 

in^ 

With values of a/h, a*, M equal to 20, .002, and 7 respectively, 
the amplification curve of Figure 34 is manifested. If subjected to 
the maximum operating pressure, buckling occurs at a time equal to 
9870 hours. If the submersible is operated at only half its operational 
depth it would not buckle until 1.85 x 10^ hours had past. 



6.5.2 SOFAR Floats 

The V7oods Hole Oceanographic Institution has employed a long 
cylinder which when released at sea sinks to a depth at which it becomes 
neutrally buoyant. The current measuring floats, first designed by 
Swallow (1963) , will follow a water mass of a constant density due to 
the neutrally compressible, buoyant properties. Placed at a depth of 
approximately 5000 feet, (the axis of the SOFAR Channel, a horizontal 
channel, in which the convex shape of the sound velocity-depth profile 
allows sound to be transmitted with minimal attenuation over thousands 
of miles) , the floats emit an acoustic tone at a given frequency. 

The float is then spatially fixed using acoustic triangulation from 
land based listening stations. The floats have been successful in 
tracking eddies in the WIIOI MODE program. Unfortunately, after two to 
three months, these floats have been found to sink an additional 230 to 
560 feet, with accelerating sinking velocity, (an apparent creep buckling 
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phenomenon) . Mavor (1974) computed a linear strain rate of 0.9 yin/in/hr 
from the change of displacement and an assumed constant in situ seawater 



density gradient. 




(6.3) 



where 



p = density 



Z = depth 



The floats are standard ALCOA extruded 6061-T6 aluminum 
cylindrical stock with a mean radius of 5.625 inches, a a/h ratio of 
7.5, an approximate out-of- roundness ratio of 0.005 and a length of 
17 feet. They operate at a mean depth of 5000 feet (2200 psi) . The 
6061-T6 aluminum, follows a cubic creep law as experimentally found by 
the author (Appendix C) 



Figure 35 displays the predicted creep response with a/h, 
a*, M equal to 7.5, 0.005, and 3, respectively. 

Since the problem is concerned with only the first three 
months of the lifetime, the perturbation transverse displacement can be 
considered linear. From Figure 35 




(6.4a) 



where 



K = 4 X 10^^ (psi) ^hr. 



(6.4b) 



w* = 1.13795 yin/hr. 



(6.5a) 



From equation (4.1), the x^erturbation hoop strain rate is 



e' = = 5.3947 x lo“^ in/in/hr. 

=2 



(6.5b) 
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From equation (4.7b), the dominant hoop strain rate is 



e, = - . (6.5c) 

0 a 

# 

Computing w from equation (4.8) yields a dominant hoop 
strain rate of 



e. = 9.45 X 10 ^ in/in/hr. (6.5d) 

U 

Hence, the total predicted linear strain rate is 

• • 

^0 ” ^0 ^*0 ” yin/in/hr, (6.5e) 

ten percent greater than the estimate of Mavor. 

This calculation shows that, in fact, creep deformation is 
most probably the cause of the unexpected sinking of the floats. 
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7. SUMMARY AND CONCLUSIONS 



The theoretical investigation which has been undertaken into 
the creep instability of the infinitely long cylindrical and complete 
spherical thin shells has satisfied the original thesis objective of 
developing a relatively simple analytical method to predict the creep 
deformation process of these shells. It was found that in both shells 
the critical buckling mode is the second elliptical mode. The roles of 
shell radius to thickness ratio and initial out-of-roundness have been 
explored. The relative superiority of the spherical shell over the 
cylindrical shell has been discussed. 

The investigation demonstrated that the single most important 
factor in creep design is designing the structure in a manner which 
maintains the stress field at the lowest possible level, thus extending 
the structure's lifetime. This one factor endorses the use of high 
strength titanium and the spherical shell in creep design. 

In conclusion, while creep buckling can be a problem in long 
^erm ocean structures, the nature of the time history of imperfection 
growth, almost negligible until the last ten percent of the lifetime, 
allows the ocean engineer to develop a safe design resistant to creep 
buckling with little difficulty. 
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8. REC0I4MENDATI0NS 



Two essential advances in creep design are required in the 
future. First, the establishment of accurate creep constants at 
ambient temperatures for common ocean engineering materials. Second, 
a more advanced theoretical investigation which encompasses the 
effects of boundary conditions of such shells as finite length cylinders, 
and partial spherical shells (end caps) . 
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w*/B _ = y' P^_ (coscf)) (cosmG) . 
m2 2 

m=0 
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- 65 - 



TITLES OF FIGURES (cont.) 



Figure 14 
Figure 15 
Figure 16 



Figure 17 



Figure 18 



Figure 19 



Figure 20 



Figure 21 



Figure 22 



Figure 23 



Figure 24 
Figure 25 
Figure 26 



Comparison of Hoff solution to herein solution for 
various m (a* = .005, a/h = 10). 

Comparison of time histories of Hoff solution to 
herein solution (m = 3, a/h =10, a* = .005) . 

Non-dimensional critical time for creep exponent 
equal to 3 for various a* and a/h for an infinitely 
long cylinder. 

Non-dimensional critical time for creep exponent 
equal to 5 for various a* and a/h for an infinitely 
long cylinder. 

Non-dimensional critical time for creep exponent 
equal to 7 for various a* and a/h for an infinitely 
long cylinder. 

Non-dimensional critical time with various initiation 
points of time segments of constant effective stress 
(i|^) (m = 3, a/h = 10, a* = .005) . 

Critical time for various initial elastic loadings, 
expressed as a percentage of yield stress, for cylinde 
of a* = .001, .01, a/h = 10, 50; and constructed of 
Titanium 621/. 08 and Aluminum 6061-T6. 

Critical time for various a/h ratios at sixty percent 
of yield stress for Titanium 621/. 08 and Aluminum 
6061-T6 cylinders (a* = .001). 

Critical time for various a* at fifty percent yield 
for Titanium 621/. 08 and Aluminum 6061-T6 cylinders 
(a/h = 10) . 

Critical time at equal external pressure for various 
a/h ratios for Titanium 626/. 08 and Aluminum 6061-T6 
(a* = .001) (P^ = 2000 psl) . 

Non-dimensional critical time for creep exponent equal 
to 3 for various R* and R/II for a complete sphere. 

Non-dimensional critical time for creep exponent equal 
to 5 for various R* and R/II for a complete sphere. 

Non-dimensional critical time for creep exponent equal 
to 7 for various R* and R/H for a complete sphere. 



TITLES OF FIGURES (cont.) 



Figure 27 



Figure 28 



Figure 29 



Figure 30 



Figure 31 



Figure 32 



Figure 33 



Figure 34 
Figure 35 



Non-dimensional critical time with various initiation 
points of time segments of constant effective stress 
(ip) (m = 3, R/H = 10, R* = .005). 

Critical time for various initial elastic loadings, 
expressed as a percentage of yield stress, for spheres 
of R* = .001, .01, R/H = 10, 50; and constructed of 
Titanium 621/. 08 and Aluminum 6061-T6. 

Critical time for various R/H ratios at sixty percent 
of yield stress for Titanium 621/. 08 and Aluminum 
6061-T6 spheres (R* = .001) . 

Critical time for various R* at fifty percent yield 
for Titanium 621/. 08 and Aluminum 6061-T6 spheres 
(R/H = 10) . 

Ratio of critical time of the complete sphere to 
critical time of the cylinder for various a/h ratios 
(m = 3) . 

Ratio of critical time of the complete sphere to 
critical time of the cylinder for various a/h ratios 
(m = 5) . ■ . 

Ratio of critical time of the complete sphere to 
critical time of the cylinder for various a/h ratios 
(m = 7) . 

Creep response of the titanium spherical hull of 
DSRV ALVIN. 

Creep response of the aluminum cylindrical SOFAR 
float. 
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APPENDIX A - COMPUTER UTILIZATION 



(i) Computer Program A - Perturbation contour map of the complete 

sphere . 

The following FORTRAN program was developed to compute . 

w*/^o for a discrete number of points, (5^ by 5^ grid) in a given 
2m 

quadrant. One needs only to input the initial values of (j) and 0 
of a .given quadrant to receive the x-y position and value of 
of the before mentioned grid. 



X = sin(j)cos0 


(1) 


y = coscj) 


(2) 


. 2 




- ^ P 2 (coscj)) cosm0 


(3) 



2m m=0 



PERTURBATION CONTOUR M7\P OF THE COMPLETE SPHERE 
DIMENSION T (19) , P(19), W(19), E(19), F(19), P2(3) 
DO 51 M=l, 4 
READ(8,200) U,V 
200 FORMAT (2F10.6) 

DO 52 N=l, 19 
A=N-1 

T (N)=.0174#5.0*A+U 
DO 53 1=1, 19 
B=I-1 

P(I)=.0174*5.0®B+V 
S=SIN(P(I) ) 

C=C0S (P (I) ) 
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4 







1 



I 



P2(l)=.5*(3.0oC*o2. 0-1.0) 



P2 (2) =3.0«S*C 

P2(3)=3.0*S**2.0 

W(I)=0.0 

DO 54 J=l, 3 

R=J-1 

D=COS (T(N)*R) 

WZ=D*P2 (J) 

54 W(I)=W(I)+WZ 

E(I)=S*COS (T(N) ) 

F(I)=C 

WRITE(5,202) E(I), P(I), W(I) 
202 FORMATC', 3F15.5) 

53 CONTINUE 
52 CONTINUE 
51 CONTINUE 
END 

Notation of Program A; 



T 



0 



P 



P2(l) 




P2 (2) 



P2(3) 



V-J 



transverse displacement/amplification function 



E 



X 
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ri 



F 



y 



u 



initial 0 value 



V 



initial (j) value 



(ii) Computer Program B - Growth of the amplification function of the 



infinitely long cylinder. 

This program was devised so that one inputs a/h, a*, M and 



ip /}p and equation (4.19) is computed. The output consists of non- 
2 1 

dimensional time (t ) , normalized amplification function (A /C ) and 

c 2 2 

a second non-dimensional time which is a better indicator of real 



INFINITELY LONG CYLINDER 

MS=3 

JS=15 

DO 1 M=l, MS 
READ(8,100) B,D,E,XM 

100 FORMAT(4F10.6) 

T=0.0 

Y=1,0 

WRITE(5,101) B,D,E,XM 

101 FORMATCl' ,5X,*A/H=* ,F10.6,’A = ' , FlO . 6 , lOX , * PHI= ' ,F10.6,10X/M 

1=‘ ,F10.6,//,15X/TIME\20X, *A1^ FN’,10X/REAL TIME',//) 

DO 2 J=1,JS 

IF (J.GT.l) GO TO 7 

F=1.0 

GO TO 8 

I 

7 y=E*x 

F=Y**(XM-1.0) 




- 103 - 



8 S=EXP (T«F/3.0) 



A=E • ( 1 . 0+6 . 0*E#B«D«EXP (T#F/3 . 0) ) 

T=3,0/F«ALOG( (A'-1.0)/(B#D^F#6.0) ) 

Tl=T/(9.0 B«*(XM+2.0)*(.866®(l.+6. B D) )• •(XM-1.)) 

WRITE (5,103) T,A,T1 

103 FORMAT (* * ,7X,F10. 3 ,17X,F10, 3 ,15X,E20. 5) 

2 X=Y 

1 CONTINUE 
END 

The number of input statements is controlled by the variable 
MS, while the number of cycles for a given set of inputs is controlled 
by the variable JS • 



Notation of Program B; 



T 



T 

C 



A 

T1 



A /C 
2 



tP 



m 



2 



Y 

XM M 



B a/h 

D a* 



E 



2 



1 



(iii) Computer Progr£un C - Growtli of the normalized spherical amplifi- 
cation function with time. 
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This program is structured exactly as the cylindrical use* 



Equation (5.1 ) is computed from the inputs R/H, R* M and ijj /ip . 

1 2 

The output received is non-dimensional time (t^) , amplification function 

(B^ /D^ ) and a second non-dimensional time which is a better indicator 
2m 2m / _ m> 



of real time 






Again, the number of input statements and number 



of cycles for a given set of input is controlled by the variables MS, 
JS respectively. 



COMPLETE SPHERE 

MS=8 

JS=50 

DO 1 M=1,MS 
READ(8,100) B,D,E,XM 

100 FORMAT (4F10. 6) 

T=0.0 

Y=1.0 

WRITE(5,101) B,D,E,XM 

101 FORMAT (' 1* ,5X, 'A/H=' ,F10.6, 'A= ’ ,F10 . 6 , lOX , ’ PHI= ’ ,F10.6,10X, 'M 

1=' ,F10.6,//,15X, 'TIME' ,20X, 'AMP FN ’ , lOX , ' REAL TIME',//) 

DO 2 J=1,JS 

IF (J.GT. 1) GO TO 7 

F=1.0 

GO TO 8 

7 Y=E«»X 
F=Y*»(XM-1.0) 

8 S=EXP(T-»F/3.0) 
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R1=4.5*B*D*F^S 
R2=112.5*(B*D«F»S)® *2.0 
G=E • »2 . 0 • (0 . 5+R1+R2 ) -0 . 5 

A=( (-'Rl/S) + ( (Rl/S)* o2.0+4.0«G®R2/So •2.0)® ®0 . 5) / (2 . 0*R2/S» *2.0) 
T=3.0/F*ALOG(A) 

Tl=T*2./(9.*B*#(XM+2. ) • (1./2.* •. 5 • ( . 5~4 . 5*B oD+112 . 5« (B«D) • *2 . ) ) 
0.5*(XM-1.) 

WRITE (5, 103) T,A,T1 

103 FORMAT ( * * , 7X ,F10 . 3 , 17X , FlO , 3 , 15X ,E20 . 5) 

2 X=Y 
1 CONTINUE 
END 



Notation of Program C; 
T 
Y 
A 

T1 

XM 

B 

D 

E 



^2m^^2m 



tP 



m 



k 

M 

R/II 

R* 
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APPENDIX B - EQUIVALENT NOTATION FROM PAST INVESTIGATIONS 



(i) Hoff^ N. J, y Jahsman^ W.E. and Nachbar^ W. "A Study of Creep 
Collapse of a Long Circular Cylindrical Shell Under Uniform 
External Pressure", J. of Aero Space Science, Vol. 20, No. 10, 



1959. 




Notation of Paper 


Equivalent Notation 


a 


a* 


X = -as 
P 


2a. (i) 


7' = (J)"^ 

3 Aa p 

1 9 9 

t = “ logd + —/X 


^ ') 


cr 2 ^ 2' Q 


2t“1 


A 


h/1 


k 

x" 


K 


n 

where P ” ^ gyration 

R = average radius 
a = shape factor 
X = weighted shape factor 


M 



(ii) Berman, 1., Chern, J.M. and Cupta , G.D. "A Parametric Study 

of Elastic-Plastic Creep Buckling of a Thin Cylindrical Shell", 
J. of Pressure Vessel Technology, August 1974. 
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Notation of Paper 



Equivalent Notation 



(iii) 



w 

a 
h 

q 

Samuelson, L.A. *'Creep Buckling o 
Non-Uniform External Loads", Int. 
1970. 

Notation of Paper 

B 

R 

The remainder of the notation is 



a* 

a 

h 

P 

s 

■ a Cylindrical Shell Under 
J. Solids Structures, Vol. 6, 

Equivalent Notation 
0 
K 

a 

;onsistent. 
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APPENDIX C - CREEP CONSTANTS OF COMMON 
OCEAN ENGINEERING MATERIALS 



(i) Titanium 621/. 08 (Ti-6AL-lTa- . 8 Mo) 

(a = 120 ksi) 

y 

Chu (1972) found that the creep data of Ti 621/. 08 at room 
temperature can be represented as 



where 



g 13.890 0.183 

4o3.55^ ^ 



( 1 ) 



e = creep strain, % 
a = stress, ksi 
t = time , hours 

If equation (1) is plotted at a given stress level, it is 
discovered that after an initial period of primary creep, the non- 
linear viscous behavior of secondary creep becomes apparent 

as example: for a = 103.55 ksi, e = .1071 yin/in/hr. 

By calculating the constant strain rates for a number of 
stress levels the following power law emerges 

= 8.3 X 10‘*°e (2) 

where a = stress, psi 

e = creep strain rate , yin/in/hr 
(ii) Titanium 6-4 (6AL-4V) 

(a = 135 ksi) 

y 

Chu has found the emperical equation for titanium 6-4 at 
room temperature, but has not as yet published it. 
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(iii) Aluminum 6Q61-T6 (wrought alloy - magnesium and silicon alloy 



solution treated and artifically aged) 

(a = 40 ksi) 

y 

The author conducted a number of tests in which cylindrical 
columns of Aluminum 6061-T6 were loaded with a uniform load in compres- 
sion by a lever system, ’Deformation was measured using two uniaxial, 
strain gages mounted on opposite sides of the column in a series con- 
figuration. The series configuration allowed bending effects to be 
internally cancelled out of the readings by the gages. 

The Aluminum 6061-T6 was found to abhor to a cubic law at 
room temperature 

= 4 X lO^^e (3) 

where 

a = stress, psi 

e = creep strain rate, yin/in/hr. 

All tests were conducted at high load levels, greater than 
65 percent of yield, at room temperature. The strain rates computed 
are 

o (psi) e (yin/in/hr) 

34670 0.8504 

33750 0.8420 

31450 0,7768 

27615 0,0569 

(iv) Stainless Steel 304 (Austentic) 

Berman (1974) expresses the secondary creep behavior of 
Stainless Steel 304 at 1050^F as 
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o 5'6789 ^ 8.035 x IO^^e 



where 



a = stress, psi 
e = creep strain, yin/in/hr. 
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